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ABSTRACT 

We show that effective chiral Lagrangians endowed with Brown-Rho scaling can be mapped 



Ph. to Landau Fermi-liquid fixed point theory in a way consistent with general constraints 



following from thermodynamics. This provides a unified scheme to treat, starting from 



^ ■ normal nuclear matter, hadronic matter under extreme conditions that is encountered in 

relativistic heavy- ion collisions and in the interior of compact stars. 
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1 Introduction 

In a recent publication [|J, we proposed an effective chiral Lagrangian whose param- 
eters (masses and coupling constants) scale in nuclear medium according to the scaling law 
described by Brown and Rho Q (referred to in what follows as "BR scaling" in short) that, 
in the mean field approximation, is to describe the ground state of nuclear matter while 
when fluctuated in various flavor directions should enable us to extrapolate to a density 
regime beyond the normal matter as well as to treat excitations above the ground state. 
The simplest such Lagrangian takes the form[| 

£ = $hp(idV-gt(p)uf)-M*(p) + hj>]il> (1) 
- m?(p)<t> 2 ] - \Fl + \m^{p)u? 

where ifj is the nucleon field, co^ the isoscalar vector field, (ft an isoscalar scalar field and the 
masses with asterisk are BR-scaling as first introduced in The scaling behavior of the 
constant g v is left arbitrary and the h is assumed not to scale although it is easy to take 
into account the density dependence if necessary. It was shown in Q that in the mean field, 
this Lagrangian - with the BR scaling suitably implemented - gives a surprisingly good 
description of the ground state with a compression modulus well within the accepted value 
200 ~ 300 MeV. 

As it stands, the Lagrangian (|l|) does not look chirally invariant. This is because 
we have dropped pion fields which play no role in the ground state of nuclear matter. In 
considering fluctuations around the ground state, they (and other pseudo-Goldstone fields 
such as kaons) should be reinstated. The chiral singlet uj field and <p field can be considered 
as auxiliary fields brought in from a Lagrangian consisting of multi-Fermion field operators 
0] via a Hubbard-Stratonovich transformation. 

The simple Lagrangian (|l|) embodies the effective field theory of QCD discussed by 

x As suggested in j^, [1, chiral in-medium Lagrangians can be brought to a form equivalent to a Walecka- 
type model. The scalar field appearing here transforms as a singlet, not as the fourth component of 0(4) of 
the linear sigma model. 
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Furnstahl et al || anchored on general considerations of chiral symmetry. As argued in Q , 
this Lagrangian should be viewed as an effective Lagrangian that results from two successive 
renormalization group "decimations" , one leading to a chiral liquid structure Q at the chiral 
symmetry scale and the other with respect to the Fermi surface [7J|. The advantage of ([!]) 
is that it can, on the one hand, be connected to Landau Fermi liquid fixed point theory 
of nuclear matter as suggested in |l|, || and, on the other hand, be extrapolated to the 
regime of hadronic matter produced under extreme conditions as encountered in relativistic 
heavy ion processes. It would, for instance, allow one, starting from the ground state of 
nuclear matter, to treat on the same footing the dilepton processes observed in CERES 
experiments as explained in || and kaon production at SIS energy and kaon condensation 
in dense matter relevant to the formation of compact stars as discussed in [K|] . 

In this note, we address an issue which was left unaddressed in namely thermo- 
dynamic consistency of the Lagrangian (|T|) treated in the mean field approximation. For 
instance, it is not obvious that the presence of the density-dependent parameters in the La- 
grangian does not spoil the self-consistency of the model, in particular, energy-momentum 
conservation in the medium and also certain relations of Fermi-liquid structure of the mat- 
ter^]. The purpose of this paper is to show that there is no inconsistency in doing a field 
theory with BR scaling masses and other parameters. 

2 Implementing BR Scaling in the Lagrangian 

In [Q] , we have treated the density-dependent masses and constants as independent of 
the fields that enter in the Lagrangian. The Euler-Lagrange equations of motion are then 
the same as for the Lagrangian wherein the masses and constants are not BR-scaling. While 
this procedure gives correct energy density, pressure and compression modulus, the energy- 
momentum conservation is not automatically assured. In fact, if one were to compute the 

pressure from the energy-density <?, one would find that it does not give | < Ta > (where 
2 We thank Bengt Friman for raising this question. 
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T^ u is the conserved energy-momentum tensor and the bra-ket means the quantity evaluated 
in the mean-field approximation as defined before) unless one drops certain terms without 
justification. This suggests that it is incorrect to take the masses and coupling constants 
independent of fields in deriving, by Noether theorem, the energy-momentum tensor. So the 
question is: how do we treat the field dependence of the BR scaling masses and constants? 

One possible solution to this problem is as follows. In Q, the density dependence of 
the Lagrangian arose as the "vacuum" expectation value]] of the scalar field % that figures 
in the QCD trace anomaly. It corresponded to the condensate of a quarkonium component 
of the scalar x with the gluonium component - which lies higher than the chiral scale - 
integrated out. It was assumed to scale in dense medium in a skyrmion-type Lagrangian 
subject to chiral symmetry. Now in the language of a chiral Lagrangian consisting of the 
nucleonic matter field ip with other massive fields integrated out, this scalar condensate 
would be some function of the "vacuum" expectation value of iptp or i[)^Qip coming from 
multi-Fermion field operators mentioned above. How these four-Fermi and higher-Fermi 
field terms can lead to BR scaling in the framework of chiral perturbation theory was 
discussed in We shall follow this strategy in this paper leaving other possibilities (such 
as dependence on the mean fields of the massive mesons) for later investigation. For this, 
it is convenient to define 

pu^ = ijrfil} (2) 
with unit fluid 4- velocity = . ^ (1, v) = , 1 = (p,j) with the baryon current density 

yl—v p 2 —j 2 

j =< ^pjip > and the baryon number density p =< ip^ifi >= J2iPi- We will take pi to be 
given by the Fermi distribution function, pi = 9(kp — \ki\) at T = 0. We should replace p 
in ([!]) by p for consistency of the model. The definition of p makes our Lagrangian Lorentz- 
invariant which will later turn out to be useful in deriving relativistic Landau formulas. 



3 By "vacuum" we mean the state of baryon number zero modified from that of true vacuum by the strong 
influence of the baryons in the system. See later for more on this point. 
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With this, the Euler-Lagrange equation of motion (EOM) for the nucleon field is 

6C d£ d£dp_ 

Sip dip dp dip 

= [i-f fl (d fl + ig*uj II -iu ll t)-M i ' + h4>}i; = (3) 



with 

^ dp 

- <^-<^-^%-^- w 

This additional term which may be related to what is referred to in many-body theory as 
"rearrangement terms" plays a crucial role in what follows. The EOM's for the bosonic 
fields are 

(#*0 M + ro* 2 )0 = hipip (5) 
W + m£V = glH^. (6) 

3 Equation of State (EOS) 

We start with the conserved canonical energy-momentum tensor constructed a la 
Noether from the Lagrangian (Q): 

T i»> = i$ 1 »d v il) + d> x <l>d v <p:-d> x u x d v uj x 

_Ip )2 _ m ^2 _ {duj) 2 + m * 2(J 2 _ ^ . u ^ r _ (?) 

We shall compute thermodynamic quantities from (]?]) using the mean field approximation 
which amounts to taking 



if) = -y=J2 a i\H^ — ~ ff-iu ] exp (iKi ■ x - i(g*u - uoZ + Ei)t) (8) 



H = C 2 h <i>il>>=ClY,Pi i ™ L 5 ( 9 ) 

i yjnf + mf 

g*u = C 2 vP = C 2 v Y.Pi (10) 

i 

= clUclY.Pi i = = (ii) 

i JHf+mf 



where is the annihilation operator of the nucleon i, with pi =< ajaj >, and £ = (£), 



= \Jk 2 + mj 2 with m£ = Af* — /i<^>, \ is the spinor and a is the Pauli matrix. We have 
defined 

C M = 4% ( 12 ) 

Cfc(/0 = -4r~, ( 13 ) 

2 = k-C 2 j+uE. (15) 

In this approximation, the energy density is 

8 = <T 00 > 

= < iifrfSPif) + ^mf(p 2 - ^m*V + tnjyfutf}) > 

= lc 2 (p 2 +r) + ^C 2 (mi-M^ 2 + ^2pi^W+<-^-l (16) 

Note that the ^-dependent terms cancel out in the comoving frame (v = 0), so that the 
resulting energy-density is identical to what one would obtain from the Lagrangian in the 
mean field with the density-dependent parameters taken as field-independent quantities as 
in H. 



Given the energy density (|i6[), the pressure can be calculated by (at T = 0) 



dE 2 d8/p 

p = -dv =p ^p- = flp - £ 

= ±C 2 v {p)p 2 -Y< p-\c 2 (p)(rnl-M*{p)) 2 

where p is the chemical potential - the first derivative of the energy density with respect 
to p in the comoving frame (v = 0): 

P = §^£\v=o = C 2 p + Ep — Sq (18) 
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with Ep = J kp + m*^ and So = (S)^ =0 . To check that this is consistent, we calculate the 
pressure from the energy-momentum tensor (0) in the mean field at T = 0: 



Pt = g < Ta >rf=o 

= \cl{p)9 2 - \cl{p){ml - M*(p)) 2 - Sop 

This agrees with (|17|). Thus our EOS conserves energy and momentum. 

4 Landau Fermi-Liquid Parameters 

The next issue we address is the connection between the mean-field theory of the chiral 
Lagrangian (|T|) and Landau's Fermi-liquid fixed point theory as formulated in |j. As 
far as we know, this connection is the only means available to implement chiral symmetry 
of QCD in dense matter based on effective field theory. For this, we shall follow closely 
Matsui's analysis of Walecka mean field model |ll|] exploiting the similarity of our model 
to the latter. 

4.1 Quasiparticle interactions 

The quasiparticle energy £j and quasiparticle interaction /jj are, respectively, given by first 
and second derivatives with respect to pf 

£i = Wi 1 fa = Wi' ( 20 ) 



A straightforward calculation gives 



(21 ) 
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and 



113 d Pj ^=° 



,2 , An 9C v , rn* L dm* L 2u dC v\2 , n 9 C v, 



+(">*L - M-n^f + C„^] + 2&^(mJ - _ M *) 



dp v * 7 dp n dp d P] y ^ ' ny h ' dp 
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^ 2 S , h dj 



-(C -^)^--^. (22) 
P Ei dpj 



with Ei = \jk 2 + m*/ 2 . Note that C v , Ch, and M* are functions of (p) = uop — u ■ j in the 
mean field approximation. In arriving at (^), we have used the observation that in the 
limit j — > 0, we have 



0, 



duo 
dpi 

d 2 u 1 dj dj 



dpidpj p 2 dpi dpj ' 

du 1 dj 

dpi p dpi ' 

^ - 1, 
dpi 

<9 2 (p) 1 dj dj 



dpidpj p dpi dpj 

and that if / is taken to be a function of the expectation value of p, then as j —* 0, we have 

df df d(p) df 



dpi d(p) dpi dp 

d 2 f _ d 2 f d(p) d(p) df d 2 (p) 



(23) 



dpidpj d(p) 2 dpi dpj d(p) dpidpj 

d 2 f Idf dj dj 
dp 2 p dp dpi dpj 

In the absence of the baryon current, j = 0, the quantities and -j^- simplify to 



(24) 



E l 
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and 



= Ei , . (26) 

dPj l + {C 2_^ )Tjipi M±< 

Writing in the standard way 

r HQ k- ■ k- 

fi = (21 + 1) J —^(^2)^.(1^1 = \ kj \ = k F ), (27) 

we see that the last term in (22) contributes to f\ and the sum of the rest at the Fermi 
surface (i.e. \kj\ = kp) to /o- So 

_ ^ykpEp 3E F 

F ° = -^ k -Tp- pk 

- 3Ef n \r 2 + AC n dCv + ^ 9m * L + n 2 \( 9Cv \ 2 + C (9R\ 

~ Tp- p[Cv + iCvp ^i + Ep-^~ + p + Cv W~ ] ( } 

Hmi - M*f[(^f + C h ^] + 2C h ^(mt - M*)±(m* L - M*) 
-2C h — (m L - M )— - C,— — (r^ - M*) - C h (m L - M )-^- 



and 



lk F E F , 3(C'l - Si)p 



1 - -ss-A = - Ef + ( c|->' (29) 



4.2 Compression modulus and F 
The compression modulus K defined by 

K S9P 5!£|^ ,30, 

comes out to be 

K - ^ + 9 ^ + 4C ^ + ^^7 + /5{( ^ } (31) 

- M*) 2 {( |^ + C^} + 2C^« - ^)|« - M*) 
2C h Qp [m L M ) ^ C h ^ Qp (m L M ) C h (m L M ) g p2 I 
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Comparing (|28|) and (pi]), we verify that our model satisfies the relativistic Landau Fermi- 
liquid formula for the compression modulus [HZ 



K = ^(l + F ). (32) 



4.3 First sound velocity 

The first sound velocity c\ in the relativistic case is defined by 



2 _ d P 
Cl = d£ 

d_fod£/p\ 
38 V dp J 
dp d 

K 



From (p2|), we have that 



(33) 



ci = v F ^(l + F ). (34) 
This is of course satisfied in our model. 
4.4 Relativistic Landau effective mass 

Baym and Chin have shown [12] that the relativistic Landau liquid satisfies the mass relation 

k F (p-) =n(l + F 1 /3). (35) 

\0£iS k=k F ,v=0 



In our model hp (ff 1 ^ ^ Q = Ep. One can see from equations (|i~8| ) and (29) that (^5|) 
is satisfied exactly in our model. 

5 Discussions 

We showed that a simple effective chiral Lagrangian with BR scaling parameters is 
thermodynamically consistent, a point which is important for studying nuclear matter under 
extreme conditions. It is clear however that this does not require that the masses appearing 
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in the Lagrangian scale according to BR scaling only. What is shown in this paper is that 
masses and coupling constants could depend on density without getting into inconsistency 
with general constraints of chiral Lagrangian field theory. This point is important for 
applying ([l]) to the density regime p ~ 3/3o appropriate for the CERES dilepton experiments 
and also kaon production at GSI where deviation from the simple BR scaling of [Q] might 
occur. 

The crucial question is really how to understand the scaling masses and constants as 
one varies temperature and density as considered in Q. If one takes the basic assumption 
of 0, | that the chiral La grangian in the mean field with BR scaling parameters corre- 
sponds to Landau's Fermi-liquid fixed point theory, then one should consider first fixing the 
Fermi momentum kp and let renormalization group flow come to the fixed points of the 
effective mass M* for the nucleon and Landau parameters T In this case, the scaling 
quantities would seem to be dependent upon A/kp, not on the fields entering into the effec- 
tive Lagrangian. This paper however shows that if one wants to approach the Fermi-liquid 
fixed point theory starting from an effective chiral Lagrangian of QCD, it is necessary to 
take into account the fact that the scaling arises from the effect of multi- Fermi interactions 
figuring in chiral Lagrangians as implied by chiral perturbation theory described in Q. 
This is probably due to the fact that we are dealing with two-stage "decimations" in the 
present problem - with the Fermi surface formed from a chiral Lagrangian as a nontopo- 
logical soliton (i.e., "chiral liquid" Q) - in contrast to condensed matter situations where 
one starts ah initio with the Fermi surface without worrying about how the Fermi surface 
is formed. Our result suggests that there will be a duality in describing processes manifest- 
ing the scaling behavior. In other words, the change of "vacuum" by density exploited in 
P] could equally be represented by a certain (possibly infinite) set of interactions among 
hadrons - e.g., four- Fermi and higher-Fermi terms in chiral Lagrangians - canonically taken 
into account in many-body theories starting from the usual matter-free vacuum. A notable 
evidence may be found in the two plausible explanations of the low-mass enhancement in 
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CERES dilepton yields in terms of scaling vector-meson masses |J and in terms of hadronic 
interactions giving rise to increased widths [jDj ]. 

How to go from one decimation to the next in hadronic physics remains an open 
problem as stressed in Q. 
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